Comprehensive genome assembly reveals genetic diversity and carcass consumption insights in critically endangered Asian king vultures.

The Asian king vulture (AKV), a vital forest scavenger, is facing globally critical endangerment. This study aimed to construct a reference genome to unveil the mechanisms underlying its scavenger abilities and to assess the genetic relatedness of the captive population in Thailand. A reference genome of a female AKV was assembled from sequencing reads obtained from both PacBio long-read and MGI short-read sequencing platforms. Comparative genomics with New World vultures (NWVs) and other birds in the Family Accipitridae revealed unique gene families in AKV associated with retroviral genome integration and feather keratin, contrasting with NWVs' genes related to olfactory reception. Expanded gene families in AKV were linked to inflammatory response, iron regulation and spermatogenesis. Positively selected genes included those associated with anti-apoptosis, immune response and muscle cell development, shedding light on adaptations for carcass consumption and high-altitude soaring. Using restriction site-associated DNA sequencing (RADseq)-based genome-wide single nucleotide polymorphisms (SNPs), genetic relatedness and inbreeding status of five captive AKVs were determined, revealing high genomic inbreeding in two females. In conclusion, the AKV reference genome was established, providing insights into its unique characteristics. Additionally, the potential of RADseq-based genome-wide SNPs for selecting AKV breeders was demonstrated.

De novo assembly and analysis of Sonneratia ovata genome and population analysis.

Mangroves are an important part of coastal and estuarine ecosystems where they serve as nurseries for marine species and prevent coastal erosion. Here we report the genome of Sonneratia ovata, which is a true mangrove that grows in estuarine environments and can tolerate moderate salt exposure. We sequenced the S. ovata genome and assembled it into chromosome-level scaffolds through the use of Hi-C. The genome is 212.3 Mb and contains 12 chromosomes that range in size from 12.2 to 23.2 Mb. Annotation identified 29,829 genes with a BUSCO completeness of 95.9%. We identified salt genes and found copy number expansion of salt genes such as ADP-ribosylation factor 1, and elongation factor 1-alpha. Population analysis identified a low level of genetic variation and a lack of population structure within S. ovata.

A de novo chromosome-scale assembly of the Lablab purpureus genome.

Introduction: Lablab (Lablab purpureus (L.) Sweet), an underutilized tropical legume crop, plays a crucial role in global food and nutritional security. To enhance our understanding of its genetic makeup towards developing elite cultivars, we sequenced and assembled a draft genome of L. purpureus accession PK2022T020 using a single tube long fragment read (stLFR) technique. Results and discussion: The preliminary assembly encompassed 367 Mb with a scaffold N50 of 4.3 Mb. To improve the contiguity of our draft genome, we employed a chromatin contact mapping (Hi-C) approach to obtain a pseudochromosome-level assembly containing 366 Mb with an N50 length of 31.1 Mb. A total of 327.4 Mb had successfully been anchored into 11 pseudomolecules, corresponding to the haploid chromosome number in lablab. Our gene prediction recovered 98.4% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Comparative analyses utilizing sequence information from single-copy orthologous genes demonstrated that L. purpureus diverged from the last common ancestor of the Phaseolus/Vigna species approximately 27.7 million years ago. A gene family expansion analysis revealed a significant expansion of genes involved in responses to biotic and abiotic stresses. Our high-quality chromosome-scale reference assembly provides an invaluable genomic resource for lablab genetic improvement and future comparative genomics studies among legume species.

Adaptation of African swine fever virus to MA-104 cells: Implications of unique genetic variations.

African swine fever virus (ASFV) is a large, double-stranded DNA virus that causes a fatal, contagious disease specifically in pigs. However, prevention and control of ASFV outbreaks have been hampered by the lack of an effective vaccine or antiviral treatment for ASFV. Although ASFV has been reported to adapt to a variety of continuous cell lines, the phenotypic and genetic changes associated with ASFV adaptation to MA-104 cells remain poorly understood. Here, we adapted ASFV field isolates to efficiently propagate through serial viral passages in MA-104 cells. The adapted ASFV strain developed a pronounced cytopathic effect and robust infection in MA-104 cells. Interestingly, the adapted variant maintained its tropism in primary porcine kidney macrophages. Whole genome analysis of the adapted virus revealed unique gene deletions in the left and right variable regions of the viral genome compared to other previously reported cell culture-adapted ASFV strains. Notably, gene duplications at the 5' and 3' ends of the viral genome were in reverse complementary alignment with their paralogs. Single point mutations in protein-coding genes and intergenic regions were also observed in the viral genome. Collectively, our results shed light on the significance of these unique genetic changes during adaptation, which facilitate the growth of ASFV in MA-104 cells.

A SNP variation in the Sucrose synthase (SoSUS) gene associated with sugar-related traits in sugarcane.

Background: Sugarcane (Saccharum spp.) is an economically significant crop for both the sugar and biofuel industries. Breeding sugarcane cultivars with high-performance agronomic traits is the most effective approach for meeting the rising demand for sugar and biofuels. Molecular markers associated with relevant agronomic traits could drastically reduce the time and resources required to develop new sugarcane varieties. Previous sugarcane candidate gene association analyses have found single nucleotide polymorphism (SNP) markers associated with sugar-related traits. This study aims to validate these associated SNP markers of six genes, including Lesion simulating disease 1 (LSD), Calreticulin (CALR), Sucrose synthase 1 (SUS1), DEAD-box ATP-dependent RNA helicase (RH), KANADI1 (KAN1), and Sodium/hydrogen exchanger 7 (NHX7), in a diverse population in 2-year and two-location evaluations. Methods: After genotyping of seven targeted SNP markers was performed by PCR Allelic Competitive Extension (PACE) SNP genotyping, the association with sugar-related traits and important cane yield component traits was determined on a set of 159 sugarcane genotypes. The marker-trait relationships were validated and identified by both t-test analysis and an association analysis based on the general linear model. Results: The mSoSUS1_SNPCh10.T/C and mSoKAN1_SNPCh7.T/C markers that were designed from the SUS1 and KAN1 genes, respectively, showed significant associations with different amounts of sugar-related traits and yield components. The mSoSUS1_SNPCh10.T/C marker was found to have more significant association with sugar-related traits, including pol, CCS, brix, fiber and sugar yield, with p values of 6.08 × 10-6 to 4.35 × 10-2, as well as some cane yield component traits with p values of 1.61 × 10-4 to 3.35 × 10-2. The significant association is consistent across four environments. Conclusion: Sucrose synthase (SUS) is considered a crucial enzyme involved in sucrose metabolism. This marker is a high potential functional marker that may be used in sugarcane breeding programs to select superior sugarcane with good fiber and high sugar contents.

Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance.

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

The complete chloroplast genome sequence and phylogenetic analysis of Heritiera fomes Buch.-Ham. (Malvales: Sterculiaceae).

Heritiera fomes Buch.-Ham. (1800) is a species of mangrove in the family Malvaceae, widely distributed in the Indo-Pacific and listed as 'endangered' (EN) on the International Union for Conservation of Nature's (IUCN) red list. We reported the complete chloroplast genome sequence of H. fomes. The genome was 168,521 bp in length and included two inverted repeats (IRs) of 34,496 bp, separated by a large single-copy (LSC) region of 88,604 bp and a small single-copy (SSC) region of 10,925 bp, respectively. The genome contained 87 protein-coding genes (PCGs), 8 rRNA genes, and 37 tRNA genes. The maximum-likelihood (ML) phylogenetic tree suggested that H. fomes is closely related to Heritiera angustata and Heritiera parvifolia with relatively high support bootstrap values of 86% and 100% with other species (Heritiera littoralis and Heritiera javanica), suggesting a relatively close genetic relationship between the five Heritiera plants. The chloroplast genome sequence provided a useful resource for conservation genetics studies of H. fomes and for phylogenetic studies of Heritiera.

Resequencing of durian genomes reveals large genetic variations among different cultivars.

Durian (Durio zibethinus), which yields the fruit known as the "King of Fruits," is an important economic crop in Southeast Asia. Several durian cultivars have been developed in this region. In this study, we resequenced the genomes of three popular durian cultivars in Thailand, including Kradumthong (KD), Monthong (MT), and Puangmanee (PM) to investigate genetic diversities of cultivated durians. KD, MT, and PM genome assemblies were 832.7, 762.6, and 821.6 Mb, and their annotations covered 95.7, 92.4, and 92.7% of the embryophyta core proteins, respectively. We constructed the draft durian pangenome and analyzed comparative genomes with related species in Malvales. Long terminal repeat (LTR) sequences and protein families in durian genomes had slower evolution rates than that in cotton genomes. However, protein families with transcriptional regulation function and protein phosphorylation function involved in abiotic and biotic stress responses appeared to evolve faster in durians. The analyses of phylogenetic relationships, copy number variations (CNVs), and presence/absence variations (PAVs) suggested that the genome evolution of Thai durians was different from that of the Malaysian durian, Musang King (MK). Among the three newly sequenced genomes, the PAV and CNV profiles of disease resistance genes and the expressions of methylesterase inhibitor domain containing genes involved in flowering and fruit maturation in MT were different from those in KD and PM. These genome assemblies and their analyses provide valuable resources to gain a better understanding of the genetic diversity of cultivated durians, which may be useful for the future development of new durian cultivars.

Assessment of the Genetic Diversity and Population Structure of Rhizophora mucronata along Coastal Areas in Thailand.

Unique and biodiverse, mangrove ecosystems provide humans with benefits and contribute to coastal protection. Rhizophora mucronata, a member of the Rhizophoraceae family, is prevalent in the mangrove forests of Thailand. R. mucronata's population structure and genetic diversity have received scant attention. Here, we sequenced the entire genome of R. mucronata using 10× Genomics technology and obtained an assembly size of 219 Mb with the N50 length of 542,540 bases. Using 2857 single nucleotide polymorphism (SNP) markers, this study investigated the genetic diversity and population structure of 80 R. mucronata accessions obtained from the mangrove forests in Thailand. The genetic diversity of R. mucronata was moderate (I = 0.573, Ho = 0.619, He = 0.391). Two subpopulations were observed and confirmed from both population structure and principal component analysis (PCA). Analysis of molecular variance (AMOVA) showed that there was more variation within populations than between them. Mean pairwise genetic differentiation (FST = 0.09) showed that there was not much genetic difference between populations. Intriguingly, the predominant clustering pattern in the R. mucronata population did not correspond to the Gulf of Thailand and the Andaman Sea, which are separated by the Malay Peninsula. Several factors could have influenced the R. mucronata genetic pattern, such as hybridization and anthropogenic factors. This research will provide important information for the future conservation and management of R. mucronata in Thailand.

Genotyping by Sequencing (GBS) for Genome-Wide SNP Identification in Plants.

Marker-assisted selection has played a pivotal role in developing several elite varieties in the past two decades. Molecular markers employed in plant breeding programs have recently shifted from microsatellites or simple sequence repeats (SSRs) to single nucleotide polymorphisms (SNPs) due to the ubiquity of SNP markers in the genome and the availability of various high-throughput SNP genotyping platforms. Rapid advances in sequencing technologies and the reduction in sequencing cost have facilitated SNP discovery in several plant species including non-model organisms with little or no genomic resources. Despite the lower cost of sequencing, genome complexity reduction approaches are still useful for SNP identification because many applications do not require every base of the genome to be sequenced. Genotyping-by-sequencing (GBS) is a quick and affordable reduced representation method that can simultaneously identify and genotype a large number of SNPs that has been successfully applied to a wide range of plant species. This chapter describes a robust two-enzyme GBS method for SNP discovery and genotyping that has been verified in non-model plant species.

A draft chromosome-scale genome assembly of a commercial sugarcane.

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum, S. spontaneum, and several other Saccharum species, resulting in an auto-allopolyploid with 8-12 copies of each chromosome. The current genome assembly gold standard is to generate a long read assembly followed by chromatin conformation capture sequencing to scaffold. We used the PacBio RSII and chromatin conformation capture sequencing to sequence and assemble the genome of a South East Asian commercial sugarcane cultivar, known as Khon Kaen 3. The Khon Kaen 3 genome assembled into 104,477 contigs totalling 7 Gb, which scaffolded into 56 pseudochromosomes containing 5.2 Gb of sequence. Genome annotation produced 242,406 genes from 30,927 orthogroups. Aligning the Khon Kaen 3 genome sequence to S. officinarum and S. spontaneum revealed a high level of apparent recombination, indicating a chimeric assembly. This assembly error is explained by high nucleotide identity between S. officinarum and S. spontaneum, where 91.8% of S. spontaneum aligns to S. officinarum at 94% identity. Thus, the subgenomes of commercial sugarcane are so similar that using short reads to correct long PacBio reads produced chimeric long reads. Future attempts to sequence sugarcane must take this information into account.

The complete chloroplast genome of Sonneratia griffithii Kurz (Lythraceae).

Sonneratia griffithii Kurz is a critically endangered mangrove species that can be found along the western coast of Thailand. In this study, we reported the complete chloroplast genome of S. griffithii. The chloroplast genome is 152,730 bp, consisting of one large single-copy (LSC) region, one small single-copy (SSC) region and a pair of inverted repeats (IRs). The LSC, SSC, and IR lengths are 87,226, 17,764, and 23,870 bp, respectively. The genome contains 113 unique genes, including 79 protein-coding, 30 tRNA, and 4 rRNA genes. The GC content of the chloroplast genome is 37.31%. The phylogenetic analysis based on 76 protein-coding genes showed a monophyletic group of S. griffithii and other Sonneratia species.

Genome assembly of the Pendlebury's roundleaf bat, Hipposideros pendleburyi, revealed the expansion of Tc1/Mariner DNA transposons in Rhinolophoidea.

Bats (Chiroptera) constitute the second largest order of mammals and have several distinctive features, such as true self-powered flight and strong immunity. The Pendlebury's roundleaf bat, Hipposideros pendleburyi, is endemic to Thailand and listed as a vulnerable species. We employed the 10× Genomics linked-read technology to obtain a genome assembly of H. pendleburyi. The assembly size was 2.17 Gb with a scaffold N50 length of 15,398,518 bases. Our phylogenetic analysis placed H. pendleburyi within the rhinolophoid clade of the suborder Yinpterochiroptera. A synteny analysis showed that H. pendleburyi shared conserved chromosome segments (up to 105 Mb) with Rhinolophus ferrumequinum and Phyllostomus discolor albeit having different chromosome numbers and belonging different families. We found positive selection signals in genes involved in inflammation, spermatogenesis and Wnt signalling. The analyses of transposable elements suggested the contraction of short interspersed nuclear elements (SINEs) and the accumulation of young mariner DNA transposons in the analysed hipposiderids. Distinct mariners were likely horizontally transferred to hipposiderid genomes over the evolution of this family. The lineage-specific profiles of SINEs and mariners might involve in the evolution of hipposiderids and be associated with the phylogenetic separations of these bats from other bat families.

Assessment of the Genetic Diversity and Population Structure of Rhizophora apiculata Blume (Rhizophoraceae) in Thailand.

Rhizophora apiculata is one of the most widespread and economically important mangrove trees in the Indo-West Pacific region. Knowledge of the genetic variation of R. apiculata in Thailand is limited. Here, we generated a whole-genome sequence of R. apiculata using the 10× Genomics technology. R. apiculata genome assembly was 230.47 Mb. Based on its genome, 2640 loci of high-quality biallelic SNPs were identified from 82 R. apiculata accessions collected from 17 natural mangrove forests in Thailand to assess the genetic diversity and population structure among them. A moderate level of genetic diversity of R. apiculata was observed. The average observed heterozygosity (Ho = 0.48) was higher than the average expected heterozygosity (He = 0.36). Two subpopulations were observed and confirmed from three approaches: population structure, PCA, and phylogenetic analyses. They corresponded to the Gulf of Thailand and the Andaman Sea separated by the Malay Peninsula. AMOVA analyses indicated that genetic variation was attributable to 76.22% within populations and 23.78% among populations. A high level of genetic differentiation between the two subpopulations (FST = 0.24, p < 0.001) was observed. This study evaluated the genetic diversity and population structure of R. apiculata, providing useful information for sustainable mangrove management in Thailand.

A Chromosome-Scale Genome Assembly of Mitragyna speciosa (Kratom) and the Assessment of Its Genetic Diversity in Thailand.

Mitragyna speciosa (Kratom) is a tropical narcotic plant native to Southeast Asia with unique pharmacological properties. Here, we report the first chromosome-scale assembly of the M. speciosa genome. We employed PacBio sequencing to obtain a preliminary assembly, which was subsequently scaffolded using the chromatin contact mapping technique (Hi-C) into 22 pseudomolecules. The final assembly was 692 Mb with a scaffold N50 of 26 Mb. We annotated a total of 39,708 protein-coding genes, and our gene predictions recovered 98.4% of the highly conserved orthologs based on the BUSCO analysis. The phylogenetic analysis revealed that M. speciosa diverged from the last common ancestors of Coffea arabica and Coffea canephora approximately 47.6 million years ago. Our analysis of the sequence divergence at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a genome-wide duplication in M. speciosa, agreeing with the report that members of the genus Mitragyna are tetraploid. The STRUCTURE and principal component analyses demonstrated that the 85 M. speciosa accessions included in this study were an admixture of two subpopulations. The availability of our high-quality chromosome-level genome assembly and the transcriptomic resources will be useful for future studies on the alkaloid biosynthesis pathway, as well as comparative phylogenetic studies in Mitragyna and related species.

The complete chloroplast genome sequence of Intsia bijuga (Colebr.) Kuntze (Fabaceae: Detaroideae: Afzelieae).

Intsia bijuga (Colebr.) Kuntze. (1891) is a threatened mangrove species, belonging to the Fabaceae family and is native to the western Pacific coast and Southeast Asia. Here, we applied short-read Illumina technology to sequence and assemble its chloroplast genome. The complete chloroplast genome is 158,363 bp in length, composed of one large single-copy (LSC) region of 87,489 bp, one small single-copy (SSC) region of 19,438 bp, and a pair of inverted repeats (IRs) of 25,719 bp. A total of 129 unique genes were annotated, comprising 84 protein-coding genes, eight rRNA genes, and 37 tRNA genes. Our phylogenetic analysis showed the placement of I. bijuga (OL699920.1) with Afzelia species within Fabaceae family.

The complete mitochondrial genome sequence of the karst-dwelling crab, Terrapotamon thungwa (Crustacea: Brachyura: Potamidae).

Terrapotamon thungwa is a new species of terrestrial long-legged crab discovered in a karst landscape of southern Thailand. Here, we report the first complete mitochondrial genome of this crab species. The mitochondrial genome size is 16,156 base-pairs (bp), including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA), and two ribosomal RNA (rRNA) genes. The AT and GC content of the mitochondrial genome sequence is 73.2% and 26.8%, respectively. Phylogenetic analysis with 26 crustacean species, based on 13 mitochondrial conserve genes, showed that T. thungwa was closely related to other freshwater crab species in the family Potamidae.

Chromosome-level genome assembly of Indian mangrove (Ceriops tagal) revealed a genome-wide duplication event predating the divergence of Rhizophoraceae mangrove species.

Mangrove ecosystems are unique, highly diverse, provide benefits to humans, and aid in coastal protection. The Indian mangrove, or spurred mangrove, [Ceriops tagal (Perr.) C. B. Rob.] is a member of the Rhizophoraceae family and is commonly found along the intertidal zones in tropical regions in Southeast Asia, southern Asia, and Africa. Here, we present the first high-quality reference genome assembly of the Ceriops species. A preliminary draft assembly, generated from the 10× Genomics linked-read library, was scaffolded using the proximity ligation chromatin contact mapping technique (Hi-C) to obtain a chromosome-scale assembly of 231,919,005 bases with an N50 length of 11,408,429 bases. The benchmarking universal single-copy orthologs (BUSCO) analysis revealed that C. tagal gene predictions recovered 95.8% of the highly conserved orthologs. Phylogenetic analyses suggested that C. tagal diverged from the last common ancestor of flat-leaf spurred mangrove [C. decandra (Griff.) Ding Hou] and C. zippeliana Blume ∼10.4 million yr ago (MYA), and the last common ancestor of genera Ceriops, Kandelia, and Rhizophora diverged from that of genus Bruguiera ∼49.4 MYA. In addition, our analysis of the transversion rate at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a recent whole-genome duplication in C. tagal. The STRUCTURE and principal component analyses illustrated that C. tagal individuals investigated in this study were the admixture of two subpopulations, the genetic background of which was influenced primarily by location. The availability of genomic and transcriptomic resources and biodiversity data reported in this work will be useful for future studies that may shed light on adaptive evolutions of mangrove species.

Assembly of a hybrid mangrove, Bruguiera hainesii, and its two ancestral contributors, Bruguiera cylindrica and Bruguiera gymnorhiza.

Mangroves are plants that live in tropical and subtropical coastal regions of the world, they are adapted to high salt environments and cyclic tidal flooding. Mangroves play important ecological roles, including acting as breeding grounds for many fish species and to prevent coastal erosion. The genomes of three mangrove species, Bruguiera gymnorhiza, Bruguiera cylindrica, and a hybrid of the two, Bruguiera hainesii were sequenced, assembled and annotated. The two progenitor species, B. gymnorhiza and B. cylindrica, were found to be highly similar to each other and sufficiently similar to B. parviflora to allow it to be used for reference based scaffolding to generate chromosome level scaffolds. The two subgenomes of B. hainesii were independently assembled and scaffolded. Analysis of B. hainesii confirms that it is a hybrid and the hybridisation event was estimated at 2.4 to 3.5 million years ago using a Bayesian Relaxed Molecular Clock approach.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests.

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.